Two-dimensional zero-gap electronic states at a magnetic field 



S.A. Ktitorov and Yu.V. Petrov 
Ioffe Physico-Technical Institute of the Russian Academy of Sciences 

This work was firstly published in 1986 [? ]. No real two-dimensional object with the zero-gap 
quasi-relativistic spectrum was known in that time. Such an object is well known now: this is 
graphene. That is why we decided to present it again as a e-print in a slightly modified form. 
A density of the two-dimensional zero-gap electronic states at the quantizing magnetic field in 
, the presence the Gaussian random potential has been calculated. The problem is reduced to zero- 

. dimensional spinor field theory using the holomorphic supersymmetric representation. The calcu- 

' lated density of states in the case of the mass perturbation has a delta function peak in the Dirac 

| point. This peak smears due to the potential perturbation. 

a : 

I I. INTRODUCTION 

A supersymmetry formalism for calculation of the two-dimensional Schroedinger electronic density of states (EDOS) 
in the presence of the random potential was developed in [l|. The problem was formulated in terms of the path 
', integrals. The supersymmetric holomorphic representation [H was used to project the state onto the lower 
Landau level. The problem was reduced to the zero-dimensional field theory after averaging over the random potential 
distribution. The EDOS was calculated exactly. This result was in accord with one obtained in [|[ for the case of the 
^ ' Gaussian white noise distribution. 

This approach is generalized in this work to the case of the zero-gap semiconductor with the quasirelativistic 
spectrum described by the Dirac equation in the vicinity of the bands crossing point, zero-gap semiconductors of the 
C . ; first kind (ZGSC-I) 0. 

The random background field is created by the distribution short-range impurities of two kinds: (i) impurities 
shifting the conduction and valence bands edges synchronously, i. e. shifting the chemical potential (potential 
impurities), and (ii) impurities shifting the bands edges in the opposite directions, i. e. modulating the gap width 
(chemical impurities). Such situation can be realized in IV- VI semiconductors solid solutions, where the band gap 
can be made vanishing varying temperature, pressure and composition. The spatial fluctuations of the composition 
play a role of the chemical perturbation, while the random violations of stoichiometry play a role of the potential 
perturbation. Small correlation radius of the random potential stems from the extremely strong lattice polarization in 
the IV- VI semiconductors: the Coulomb tail e 2 / (er) can be neglected, when e~400. In the two-dimensional systems 
like graphene a potentials shifting the chemical potential and modulating the mass (gap) appear because the Kohn- 
luttenger matrix elements of the short-range impurity potential are not equal generically for the upper and lower 
bands. 

A character of the EDOS is in accord with [l[ in the case of Coulomb impurities. In the case of the chemical 
impurities there is a delta function singularity in the centre of the band. Simultaneous effect of the chemical and 
Coulomb perturbations smears this singularity, when chemical perturbation dominates. When intensities of these 
perturbations are equal, EDOS is similar to one obtained in 
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II. HOLOMORPHIC REPRESENTATION FOR ZGSC-I 



■ Electronic states in ZGSC-I with impurities at the quantizing magnetic field are described by the two-band Hamil- 
^ | tonian Q: 

_ _ i 

H=H + V(r), (1) 



Ho =sa (p-^A) (2) 

where a is the Dirac matrix, s is the Fermi velocity, V (r) is the potential of impurities. We choose the axial gauge 
for the magnetic field potential: A = |H x r, H || e z . The spectrum of the operator ([2]) is well known Let us 
construct the relativistic holomorphic representation. It is carried out by functions of the following form: 



*j ( x : v) = ex P h«] u i M > 



(3) 
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where ? = (x + iy) /21h, ? = (x — iy) /2Zjy, ^ 2 = eH/hc is the magnetic length, j is the spinor index, Uj (?) is 
a holomorphic function of the variable ? in the Fock-Bargmann space We can make use of the close relation 
between the holomorphic representation and the Glauber coherent state representation The wave functions of 
the stationary states with the energy E, spin component a z = 1 and momentum k z has the following form in the 
relativistic holomorphic representation [9(: 



V2 



hsy/n/ (l H ) \n 
\n,P) 




(4) 



Here n is the Landau band number, Re/3 and Im/3 are the circular motion coordinates. It is seen from (??) that the 
small components of the bispinor vanish at n = 0. The Dirac bispinor (??) reduces into the spinor 



*o,<7,=i,fe* 0,2/) = -j= 



|0,/3) 



where |0, (3) is an eigenfunction of the operator b : 

b\O,0) =j9|0,j9) ; |0,jQT> = (V§»rZjr) * exp [-?? + v 7 ^? - |/3| 2 /2 



(•5) 



(6) 
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The function |0,/3) is a generating one for the states with a definite angular momentum z— component quantum 



number m. Let us expand exp 



exp 



0, (3) to powers of (3 : 

^max 

|0, 13) = exp (-??) B„s m = exp (-??) u (?) = 95 (a;, y) , 



(8) 



771 = 



where Z? m = — !=— /3 m /m!, m max = L 2 /2ttI%. Thus, one can see from ([5]) and (JSJ) a correspondence existing between 
the coherent and relativistic holomorphic representations. The Landau zero-band set membership is determined by 
the Cauchy-Riemann conditions for u (?) : 



«9u(?)/d? = 0, 9u(f)/9? = 0. 



(9) 



While quantizing we replace the variables ?, ? with the ladder operators ?, ? = 9/9? acting in the Fock space. Their 
commutaters read: 



(10) 



III. FUNCTIONAL INTEGRAL REPRESENTATION FOR GREEN'S FUNCTIONS 



Our goal in this section is to derive a general expression for EDOS for the quasirelativistic states of the lower 
Landau band of a two-dimensional system at the random impurity field. We follow to the approach developed in [l[ 
generalizing it to the quasirelativistic system. A general expression for EDOS reads: 



p(E) = — Im Tr [ dxdy (G (E + i0; r = r')) 
ttS J 



(11) 



where G (E + iO; r, r') is a one-particle Green's function; trace is taken on the spinor indices, S is the surface area; 
the angle brackets indicate averaging on the random potential configurations. The Green function can be presented 
by the functional integral: 



G a j,a'ji (E + i0; r, r') = -iZ 1 / DpDtp*ip aj (r) ip* a , ■, (r') exp 



i J dxdyip* aj (r) (e-H + U 



<p aj (r) 



(12) 
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D(pDip* exp 



dxdyip* aj (r) (e - H + iOj ip aj (r) 



where H. is determined by |T]). The impurity potential is given by the formulae: 

V(r) = V 1 (r)+a 1 V 2 (v), (V 4 (r)> = 0, (VJ (r) V,- (r')) = Ai%5 (r - r') 

Here V\ (r) is a screened Coulomb potential, Vi (r) is a mass (gap) modulating perturbation; o\ ■ 
matrix. 

Let us define a supervector belonging to the zero Landau level subspace: 



1 

1 



(13) 

(14) 
is the Pauli 



» i>3 = exp (-??) Wj (0 , $ = , 



where Uj (?) is the Grassmann variable. Then the averaged Green function reads in terms of the supervectors: 

Gij = -i J (?) u* (?) exp (-??) exp (i^^) , 

where the effective action reads 

Shot =E d?d?exp (-??/2) + / d?d?exp (-??) [Ai (¥$) + A 2 ($cti$) 



(15) 



(16) 



(17) 



Here Ai and A2 are intensities of the impurity correlators (|14p. It is seen from (|17p that varying of the free action 
with respect to $ gives the equation for the Landau gapless zero-band; in the two-dimensional case it is a zero-energy 
state 



Ho<P = 0. 

Following [H, we introduce the holomorphic superficlds: 



Xi = Ui (?) + A= S i V i (?) 1 Xi = < (?) + -^«< (?) 



(18) 



(19) 



where 9i are Grassmannian algebra generators. The effective action presented in terms of the superfields Xi is 
supersymmetric similarly to the Schroedinger case [l[, i. e. it is invariant under the magnetic translation group 
(translation x gauge transformation) for the Dirac electron Q and the superspace rotation group: 



S = i2irl 2 H J dOdOdsds exp 
2nl 2 H / dOdOdsdsexp 



~2 I?? 



~2 I?? 



Ai (XtXif + A21 (x l o-iX») 2 



(20) 



Xi 0) = X» ( z - a,0 - w) exp 
<5z = ZU0, <50 = ojz. 



6ZJ 



(21) 



Substituting ([20| into (TT2|) and calculating the Green function perturbatively we can see that the superspace 
Gaussian integral equals unity in all orders, while the corrections to the bare Green function are nothing but the 
zero-dimensional spinor A^ 4 field theory symmetry coefficients. The total Green function can be written as a ratio of 
ordinary Rienannian integrals; the resulting formula differs from the derived in [l[ one only by the spinor structure 
presence. 
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IV. EDOS CALCULATION 



A general formula for calculation of the EDOS for ZGSC-I at the quantizing magnetic field in the presence of the 
random imurity field reads: 



P{E) 



1 d f°° r 2 

- — - Im — log / dxdy exp iae (x + y) — a\\ (x + y) — olX^ (x — y) 2 
2-Kb oe J L 



(22) 



where a — 2itl 2 H , x ~ \u\\ 2 , y = |z/2 1 2 , Ui are spinor components, e = E + iO. Notice that the total number of 
unperturbed states per the unit area resulting from (|22[) is twice larger, than in the one-band case /_ dEp^H (E) = 
l/2nl 2 H . We will show below by the direct calculation that a presence of the chemical impurities modulating the mass 
(gap) leads to the singularity of EDOS at E = 0. Assuming Ai << A2 we can one of the integrals: 



1 d 

P ( E ) = —5 Im jr log 
nS Oe 



y / 7r/4aAi exp (— ae 2 /4Ai) erf c ^—i^y/a/Xij + 
J dxexp (2iaex) erf c \ x-\f 01X2 + i — \/a/^2j 



(23) 



We have neglected the term Aa\\x 2 in the exponent since convergence of the integral at the upper limit is guarantied 
by the complementary error function (erf c(x) — > at x — > 00 within the domain |argz| < tt/4). The integral in 
can be calculated exactly [loj : 

_ 1 (tt/4) exp (- £ 2 /A 2 ) + exp (e 2 ) [F (e/\) F (e) (l - A 2 /2 £ 2 ) - ± (F (s) - F (e/A))] 

P[E) = V^/Al^^; ; — h 



2tt 2 / 



H 



(V^F/2) exp (- £ 2 /A 2 ) - ^_ exp (s 2 ) F (s) + [F (e)] 



2irl 2 H 



6(E) 



(24) 



where e = E ^ ^ > ^ = \/^i7^2> E (x) = exp (—a; 2 ) J* dtexpt 2 . 

Let us consider some limiting cases. If A2 = 0, i. e. chemical impurities modulating the mass (gap) are absent, we 
can simplify 



P (E) = ^2/\~ 1T 1 



exp (rj 2 ) 



4tt 2 I h 2 



^ eX p( V 2 )^F( V ) 



(25) 



where 77 = -f yjaj\\. 

Taking account of the delta-correlated chemical impurities only (Ai =0, A2 7^ 0) gives a delta function peak at 
E = on the smooth background: 



p (E) = vW2A2 



1 7T 



.He) 



exp (e 2 ) jypH 



1 



^ exp (e 2 ) F (e) 



(26) 



Notice that in the case of Ai = A 2 = w/8 we obtain EDOS similar to obtained in (??), but with the additional 
facto 2: 



P {E) = ^kJ 



w 



exp (i/ 2 ) 



TT 2 l H 



(27) 



v = E^/aj 



w. 



In the limit of the Coulomb impurity correlator low intensity at 4Ai/aTi? 2 << 4A2/Q; and assuming E — > 0,we 
obtain from 



p (E) = v / 2tt/A ] 



1 



/7rcxp (e 2 ) 



8n 2 l, 



H F 2 (e/A)+ [^exp(-e 2 /A 2 ) 



(28) 
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It is seen that the peak width is proportional to y/2\\/a. When E 2 » 4X2/ a » iXi/a, we have a resuld similar 
to obtained in (??) and (??) for the one-band semiconductor: 

p{E) = ^™ti3gm. (29) 

4A 2 

Such asymptotic was obtained in P3 by the semi-classical quantization approach. 

Narrow peak appear in the strongly irregular semiconductors with the Lorentz distribution. The effective action 
takes the form S — —Ax. This gives a generalization of the Lloyd model [ll| to the zero-gap semiconductor. The 
EDOS takes the form 

where A± = \\ ± A 2 . When A_ — > (close values of intensities), the peak width tends to zeo. 



Al + E 2 A 2 _+ E 2 



(30) 



V. DISCUSSION 



In the absence of dynamic interactio the action (|20p is invariant with respect to the supersymmetry transformations. 
It is doubly degenerate (apart from the usual Landau degeneracy in the magnetic field); E = 0. Impurities do not 
violate the supersymmetry that results in appearing of the delta-peak. Possible physical realization: surface states, 
states in the supersymmetric interface in the heterojunction of mutually inverted narrow-gap semiconductors [l3l |. 
and, now the most interesting, graphene. 
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